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O . ABSTRACT 

' We present a model independent analysis of dark matter (DM) both decoupling 

ultrarelativistic (UR) and non-relativistic (NR) based on the DM phase-space density 

■ V = pdm /&dm- We derive explicit formulas for the DM particle mass m and for 

the number of ultrarelativistic degrees of freedom at decoupling. We find that 

Q . for DM particles decoupling UR both at local thermal equilibrium (LTE) and out 

' of LTE, m turns to be at the keV scale. For example, for DM Majorana fermions 

c/2 \ decoupling at LTE the mass results m ~ 0.85 keV. For DM particles decoupling NR, 

y/m Td results in the keV scale (Td is the decoupling temperature) and the m value 
is consistent with the keV scale. In all cases, DM turns to be cold DM (CDM). Also, 
^f-} , lower and upper bounds on the DM annihilation cross-section for NR decoupling are 

■ derived. We evaluate the free-streaming (Jeans') wavelength and Jeans' mass: they 
£S| result independent of the type of DM except for the DM self-gravity dynamics. The 
CN| , free-streaming wavelength today results in the kpc range. These results are based 
0^ ■ on our theoretical analysis, astronomical observations of dwarf spheroidal satellite 

galaxies in the Milky Way and ./V-body numerical simulations. We analyze and discuss 
_ . the results on T> from analytic approximate formulas both for linear fluctuations and 

(^i the (non-linear) spherical model and from TV-body simulations results. We obtain in 

Q\ \ this way upper bounds for the DM particle mass which all result below the 100 keV 

■ range. 
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1 THE DARK MATTER PARTICLE MASS and the comoving momentum of the dark matter particles 

Pc- 

Although dark matter was noticed seventy-five years ago „ . ,, e .. n , \ 

i . Imonl U J | ln/ ,nh •+ + ■ * + i Knowing the distribution function F d (p c ), we can com- 

Zwicky 1933; Oort 1940 its nature is not yet known. , . , ., , , .' a 

' , — ~ y pute physical magnitudes as the DM velocity fluctuations 

Dark matter (DM) must be non-relativistic by the time 

of structure formation (z < 30) in order to reproduce the 



(nw = (=)(()^ a) 



and the DM energy density. For the relevant times t dur- 
ing structure formation, when the DM particles are non- 
observed small structure at ~ 2 — 3 kpc. , ,. . ,. , 

r relativistic, we have 

DM particles can decouple being ultrarelativistic (UR) 

at Td 3> m or being non-relativistic (NR) at Td < m ^ J d Pyh P P h p^ a Mp j 

where m is the mass of the dark matter particles and Td the /V/ 2 \ (+\ — /^ ph — •* C^ n ^ i 

decoupling temperature. We consider in this paper particles m 2 f d s p p h 

that decouple at or out of local thermal equilibrium (LTE). 

The DM distribution function Fd freezes out at decou- , , . , , c , , . . 

where we use the physical momentum of the dark matter 

pling. Therefore, for all times after decoupling F d coincides particleg (t) _ pc/a{t) &g integration variable The gcale 

with its expression at decoupling. Fd is a function of Td, m c . ,.\ ■ >• j i 

r r ° factor a(t) is normalized as usual, 

a{t) = j + z (t) ' a ( toda -y) = 1 < ( 2 ) 

* devega@lpthe.jussieu.fr namely, the physical momentum p p h(t) coincides today with 

f Norma.Sanchez@obspm.fr the comoving momentum p c . 
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We can relate the covariant decoupling temperature Td, 
the effective number of UR degrees of freedom at decoupling 
gd and the p hoton temperature today Tj by us i ng entrop y 
conservation dKolb fc Turnenll99d ; lBornerll2003l ; lYad f2006i: 



T d = 



where = 0.2348 meV 



(3) 



and 1 meV = 10 -3 eV. 

The DM energy density can be written as 

POM(t) = g I F d [a(t)p ph ] , (4) 

where g is the number of internal degrees of freedom of the 
DM particle, typically 1 ^ g ^ 4. 

By the time when the DM particles are non-relativistic, 
the energy density eq.(|4]) becomes 



pDM (£) 



mg Tl 



2 ty 2 a 3 (t) 



m n(t) 



(•») 



where 



Jo 



V F d {y) dy , 



v ■■■ 



(6) 



n(t) is the number of DM particles per unit volume and we 
used as integration variable 

Pph(t) _ Pc_ 

Td(t) ~ T d ' 

From eq.([5|l at t = and from the val ue observed today for 
Pom |Komatsu et alJl2009l ; lYaoll200d ). 

Pdm = Qdm Pc — 0.228 p c , 

p c = 3 M 2 Pl Hi = (2.518 meV) 4 , (7) 
and M 2 Pl = 1/[8tt G], we find the value of the DM mass: 

m = tt 2 Qdm £ = 6.986 eV -if- , (8) 
T-j g h g 1 2 

where p c is the critical density. 

Using as integration variable y [eq.JBJ], eq.JTJ) for the 
velocity fluctuations, yields 

2 r 



(v 2 )(t) 



T d 



where 



m a(t) 
V 4 F d (y) dy 



(9) 



Expressing Td in terms of the CMB temperature today ac- 
cording to eq.([3]) gives for the one-dimensional velocity dis- 
persion, 



ODM (z) 



= 0.05124 



§<V 2 )(z) 

1 + z keV 



23 1 + z T- 
7! 



y] 



n 
u 
h 



km 

s 



(10) 



(11) 



It is very useful to consider the phase-space 
density i nvariant under the unive rse expan- 
sion (Boyanovsky, dc Vega & Sanc hej 2008a; 

iHogan fc Daka nton 2000; Madscn 1990, 200j) 



2>(f) = 



n(t) 



1 



ppujt) 
3 ^3m 4 (T 3 DM (t) 



(12) 



where we consider the relevant times t during structure for- 
mation when the DM particles are non-relativistic. T>(t) is 
a constant in absence of self-gravity. In the non-relativistic 
regime T>(t) can only dec rease by collisionless phase mixing 
or se lf-gravity dynamics (|Lvnden-Bellll 19671 ; iTremaine et al.l 
Il986h . 

We derive a useful expression for the phase-space den- 
sity V from eqs.©, (fT0|) and (JT2J) with the result 



V = 



In 



2 7T 2 A 



(13) 



Observing dwarf spheroidal satellite galaxies in the Milky 
Way (dSphs) yields for the phase-space density today 
<|Wvse fc Gilmorell2007l ; iGilmore et. aLlboOTT ): 



„ 5 x 10 ; 



keV/cm : 



(km/s) 



r = (0.18 keV) 4 



(14) 



The precision of these results is about a factor 10. 

After the radiation dominated era the phase-space den- 
sity reduces by a factor that we call Z 



23(0) = ^5(z~ 3200) 

Zj 



(15) 



Recall that D(z) = Pdm/{3 v3 m 4 a DM ) [according to 
eq. (|12[) ] is independent of z for z > 3200 since density 
fluctuations wer e < 10 -3 before the matter dominated era 
l|Dodelsonll2"003T ). 

The range of values of Z (which is necessarily Z > 1) is 
analyzed in detail in sec. 12.31 below. 

We can express the phase-space density today from 
eqs.JTSJl and JTUl as 



23(0) = 



3 V3 m 4 erf 
Therefore, eas.(fT2]l. JUJ and (QUI yield, 



£l = l££^(^ 3200) 
erf Z at,.. 



(16) 



(17) 



where Pdm/o- dm (z ~ 3200) follows from eqs.gl]) and JTS 



Pdm 



(z ~ 3200) = 



DM 



3 %/3 TTi 
2 7T 2 



5 — 



(18) 



We can express m from eqs. (|14)1 - (|18[1 in terms of T> and 
observable quantities as 



Z p s = _2_7r^_ Z p^I£_ 
3 n/3 © erf ~ 3 73 c/ CT 3 | 



m = 0.2504 



4r keV . 



(19) 



(20) 



Combining this with eq.(|SJ) for m we obtain the number of 
ultrarelativistic degrees of freedom at decoupling as 



yd 



2* gi T 3 fZp s 



3s 7f2 OdM Pc 

35.96 zi g'i [h 1$ 



(21) 



If we assume that dark matter today is a self-gravitating gas 
in thermal equilibrium described by an isothermal sphere 
solution of the Lane-Emden equation, the relevant quantity 
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characterizing the dynamics is the dimensionles s variable 
jde Vega fc Sanchez! |20o3 ; I Destri fc de Vegall2007h 



Gm 2 N 2 t2 Ps 
LI A at 



(22) 



which is bound t o be r\ < 1.6 to prevent the gravitational col- 
lapse of the gas (|de Vega fc Sanchedl2002l ; iDestri fe de Vegal 
l2007h . Here V = L 3 stands for the volume occupied by the 
gas, N for the number of particles, G for Newton's constant 
and T = | m a 2 is the gas tempera ture. (The length L 
is sim ilar to the so-called King radius l)Binnev fc Tremaind 
1 19871 ). Notice however that the King radius follows from 
the singular isothermal sphere solution while L is the 
characteristic size of a stable isothermal sphere solution 
l|de Veea fc Sanchez! [20021 : iDestri fc de Vegall2007l U 

The compilation of recent photomet ric and kinematic 
data from ten Milky Way d Sphs satellites l|Wvse fc Girmorel 
120071 ; iGilmore et. al.ll2007h yields values for the one dimen- 
sional velocity dispersion a s and the radius L in the ranges 

0.5 kpc ^ L 1.8 kpc , 6.6 km/s ^ a s ^ 11.1 km/s . 

(23) 

Combining eq. (ll2p . eq.JTTJ) and (|22p yields the explicit ex- 
pression for the DM particle mass, 

4 1 Z 4 71" »r2 Z 

= 0.5279 MT* ^ (keV ) 4 . (24) 

This formula provides an expression for the DM particle 
mass independent of eq. (|19[) . We shall see below that 
eqs. (|19[) and (|24|l yield similar results. 

We investigate in the subsequent sections the cases 
where DM particles decoupled UR or NR both at LTE and 
out of LTE. We compute there m and gd explicitly in the 
different cases according to the general formulas eqs, (l20|) , 
(I2T1) and HO. 



1.1 Jeans' (free-streaming) wavelength and Jeans' 
mass 

It is very important to evaluate the Jeans' length and Jeans ' 
mass in the present context ijBornerl 120031 ; I Gilbert! 1 19681; 
iBond fc Szaiavlll983r i. The Jeans' length is analogous to the 
free-streaming wavelength. The free-streaming wavevector 
is the largest wavevector exhibiting gravitational instability 
and characteri zes the scale of suppression of the DM tra nsfer 
function T(k) l|Bovanovskv. de Vega fc Sanche"3l2008bh . 

The p hysical free-st r eaming wavelength can be ex- 
presse d as (jBornerl 120031 ; iBovanovskv. de Vega fc Sanchez! 
2008b) 



X f3 (t) = \j(t) = 



2tt 
*/.(*) 



(25) 



where kf s (t) = kj(t) is the physical free-streaming 
wavenumber given by 



[! + *(*)] 



H Qdm 



2 , _ 4 7T G PDMjt) _ 3 

(V 2 )(t) ~ 2 (V 2 )(0) 



(26) 



where we used that pDMif) = Pdm(0) (1 + z) 3 and eq.©. 



We obtain the primordial DM dispersion velocity ctdm 
from eqs. ([5j), © and (|T7) |. 



(27) 



i (V 2 )(0) = gum = (3 Mpi Hi Q DM \ — 
3 V Z ps 

This expression is valid for any kind of DM particles. Insert- 
ing eq. (|27[) into eq. (|26[l yields for the physical free-streaming 
wavelength 



2V2-K f3M 



Q 6 

DM 

16.3 



Ho 



Z p s 



1 



kpc 



(28) 



(kpc)" 



where we used 1 keV = 1.563738 10 2 

Notice that A/ s and therefore Xj turn to be indepen- 
dent of the nature of the DM particle except for the factor 
Z. 

The approximated analytic evaluations in sec. 
[ j] together with the results of JV-body simulations 

ilPeirani et alj l2006t iHoffman et ail 120071 : ILapi fc Cavalierel 
120091 ; iRomano-Diaz et all 120061 . 1200?! ; IVass et all l2009h 

indicate that for dSphs Z is in the range 

1< Z < 10000 . 

Therefore, 1 < Z$ < 21.5 and the free-streaming wave- 
length results in the range 

1 , . , > „ 1 



0.757 



kpc < X f s(z) < 16.3 



kpc 



y/l + Z " c ' " J ° v " / ' y/l + z 

These values at z = are consistent with the iV-body 
simulations reported in iGao fc Theunsl (|2007t ) and are of 
the order of the small DM structures o bserved today 
l|Wvse fc Gilmore!l2007l ; IGilmore et. al]|2007D . 
The Jeans' mass is given by 



Mj(t) = - 7T X 3 j(t) PDM (t) 



(29) 



and p rovides the smal l est un s table mass by g ravitational col- 
lapse [KoIbI&~Turne3 (|l990h : lB5rnerl ^2003 ). Inserting here 
eq.© for the DM density and eq.(|28]l for Xj(t) = X fs (t) 
yields 



Mj(z) = 192 V2 7T 4 V^dm MU H 



Z P s 



Q + z) 



0.4464 



10 7 Mq (l + z)i 



(30) 



Taking into account the Z- values range yields 

0.4464 10 3 M e < Mj(z) (1 + z)~^ < 0.4464 10 7 Mq . 

This gives masses of the order of galactic masses ~ 10 Mq 
by the beginning of the MD era z ~ 3200. In addition, the 
comoving free-streaming wavelength scale by z ~ 3200 

3200 x X fs (z ~ 3200) ~ 100 kpc , 

turns to be of the order of the galaxy sizes today. 



2 THE PHASE-SPACE DENSITY V FROM 
ANALYTIC APPROXIMATION METHODS 
AND FROM iV-BODY SIMULATIONS 

We analytically derive here formulas for the reduction factor 
Z defined by eo. (|15[) in the linear approximation and in the 
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spherical model. The results obtained (see Table I) are in 
fact upper bounds for Z. We then analyze the results on 
T> from JV-body simulations. 

2.1 Linear perturbations 

The simplest calculation of T> follows by considering linear 
perturbations around the homogeneous distribution Pdm(z) 
as 

p = p D u(z)[l + 5(z,k)] , (31) 

We have Pdm(z) = Pdm(O) (1 + z) 3 and in a matter domi- 
nated universe 

1 + 2 



5{z, k) ~ Si 



1+Zi 



(32) 



The peculiar v elocity in the MD universe behaves as 
!|DodelsoJ2003h . 



aH 5~ l/vT+~ 



(33) 



We can thus relate the phase-space density at redshift z 
T>{z) ~ p/v :i eq.(fT2J) with the phase-space density at redshift 



Zi as, 



T>{z) ~ V{ Zi ) 



1 + 2 

1 + Zi 



(34) 



Since the linear approximation is valid for |<5| -C 1, we find 
from eq. (|32p that eq. (|34p applies in the redshift range (zi, z) 
where 



1 + z > (1 + Zi) Si 



(35) 



We can apply eq. (|34p to relate T>(z) at equilibration (z = 
Zi — 3200) and T>(z) at the beginning of structure formation 
2 ~ 30, since Si ~ 10 -3 as a scale average of t he density 
fluctu ations at the end of the RD dominated era ijDodelsonl 
l2003h and eq. (|35[) is satisfied. We thus obtain from eq.i 
in the linear approximation: 



V{z ~ 3200) 
V{z ~ 30) 



1.3 x 10 a 



(36) 



Notice that eq. (|35p does not hold for z ~ 0. Therefore, in 
order to evaluate Z, we should combine the linear approxi- 
mation result ea, (|36p for 30 < z < 3200 with the results of 
iV-body simulations for < z < 30. This is done in sec. 12.41 
to obtain upper bounds for Z. 

2.2 The spherical model 

Let us now consider the spherical model where particles 
only move in the radial direc tion but where the non- 



linear evolution is exactly solved ( Fillmore & Goldrich 1984 
lBertschingeilll985l : |Peebles!ll993l : |Padmanabhanlll999l ). Th< 
proper radius of the spherical shell obeys the equation 



R : 



G M 



(37) 



where G is the gravitational constant and M the (constant) 
mass enclosed by the s hell. Eq . (1371) can be solved in cl ose 
form with the solution (|Bertschingerlll985l ; |Peeb"Ieslll993f ) 

t= 1% (0-sin0) , 
4 5? 



fl=^; (1-oob( 



2R'j 
9t? 



= G M 



R : 



2Ri y/6~i sm ( 



3t 

l + z = (l + Zi)Si 



1 — cos 9 ' 
4\ § 



sin 0) 3 



Pdm{z) 



9 (e-sin< 



(38) 



2 (1 -cos0) 3 

Here, Ri and Zi are the radius and the redshift at the initial 
time U and 8 is an auxiliary time dependent parameter. 

Choosing the initial time by equilibration with 2, > 1 
we have 8i <C 1 and we find from eqs. (|38l) . 

2R, 



2 VS~ , R(ti) 



3U 



Pi = pDM(zi) . (39) 



The spherical shell reaches its maximum radius of expan- 
sion R m = Ri/Si at 8 — n and then it turns around and 
collapses to a point at 8 = 2 n. However, well before that, 
the approximation that matter only moves radially and that 
random particle velocities are small will break down. Actu- 
ally, the DM relaxes to a virialized configuration where the 
velocity and the viria l radius follow from the virial theorem 
|Padmanabhanll 19991 ) 

- ,ir ''- W ^ - 1 " (40) 



5R„ 



Rv — — Rn 



We can now compute the initial phase-space density (at Zi) 
and the phase-space density at virialization. We get at Zi 
from eqs.JT2j and J39J 



(41) 



V, 



PDM{0) (1 + 2,) 3 



2R t 



and at virialization for 9 — 2ir from eqs. (|12p and (|40p 



3 ^3 ; 



Ri 



32 
9tt 2 



15 



5, 



(42) 



Therefore, the Z factor in the spherical model takes the 
value 

3 



z = 



Vi 9^ 



32 



5 Si 



1.29009 



Setting Si ~ 10 3 as a scale average of the density fluctua - 
tions at the end of the RD dominated era (|Dodelsodl2"003h 
yields 

->4 



4.08 x 10* 



(43) 



The spherical model approximates the evolution as a purely 
radial expansion followed by a radial collapse. Since no 
transverse motion is allowed neither mergers, the spherical 
model result for Z eq. (|43p is actually an upper bound on Z. 



2.3 The phase-space density V from N body 
simulations 

The phase-space density T>(z) is invariant under the 
universe expansion except for the self-gra vity dynam- 
ics that diminishes T>(z) in its evolution (|Lvnden-Belll 
1 19671 : iTremaine et al.1 1 19861 ). Numerical simulations show 
that T>(z) decreases sharply during phas es of violent 
mergers followed by quiescent phases |Peirani et al.l 
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Approximation used 


Upper limit on Z 


i 

Upper limit on m ~ 0.5 Z & keV 


Linear fluctuations 


~ 1.3 X 10 11 


96 keV 


Spherical Model 


3 

~ 1.29 X S i 2 ~ 4.1 x 10 4 


7.1 keV 



Table 1. Upper bounds for the Z-factor [defined by eg. 1151 1 and for the mass of the DM particle obtained 
for two different approximation methods. Notice that only the spherical model takes into account non-linear 
self-gravity effects. The mass m mildly depends on Z through the power 1/4. In any case m results in the 
keV range. 



20061; iHoffman et al.1 [20071; lLapi fe Cavalierd l2009t 
Romano-Diaz et al l 12009 . 120071 ; IVass et all 1200a) . V{z) 



decreases at these violent phases by a factor of the or- 
der > 1. (Se e fig 3 in IPeirani et all (I2006TI , fig. 1 in 
IHoffman et ail J2007T), fig 6 in lLapi fe Cavalierd (|2009h and 
fig. 5 in IVass et al.l (|2009r )). These sharp decreasings of T> 
are in agreement with the linear approximation of sec. 12.11 
as we show below. 

A succession of several violent phases happens during 
the structure formation stage (z < 30). Their cumulated ef- 
fect together with the evolution of T> for 3200 > z > 30 
produces a range of values of the Z factor which we can 
conservatively es timate on the basi s of the A^-body sim- 
ulations results (IPeirani et alj 20061; Hoffman et al. 2007 ; 



Lapi fe Cavaliere! 120091 ; iRomano-Diaz et al.l 120061 . 12007 ; 
Vass et al.ll2009l ) and the approximation results eos. (|36p and 



((43)) . This gives a range of values 1 < Z < 10000 for dSphs. 

Indeed, more accurate analysis of N body simulations 
should narrow this range for Z which depends on the type 
and size of the galaxy considered. 

The dSphs observations, for which the best observa- 
tional data are available, take mostly into account the cores 
of the structures since the dSphs have been stripped of their 
external halos. Hence, the observed values of p s /o"s may be 
higher than the space-averaged valued represented by the 
right hand side of eqs. (|15|) and (|17|) . Higher values for p B /al 
correspond to lower values for Z. 



The approximate formula eg. (| 34 1) indicates a sharp 
decrease of the phase-space density with the red- 
shift. This sharp decreasing is in qualitati ve agreement 



with th e simulations in the violent 
(2006)}; IHoffman etal 



pha 



., i 2007li; lLapi fc Cavalier. 

IRomano-Diaz et al l (|200d . |2007» : IVass et al l (| 20p" 



Peirani et al.l 
(|2Q09D; 



2.4 



Sinthetic discussion on the evaluation of Z 
and its upper bounds. 



The DM particle mass scale is set by the phase-space density 
for dSphs eq. (|14[) . Those galaxies are particularly dense and 
exhibit larger values for Ps/crf than spiral galaxies. Since the 
primordial phase-space density is an universal quantity only 
depending on cosmological parameters, the ^-factor must be 
galaxy dependent, larger for spiral galaxies than for dSphs. 

We want to stress that the values of the relevant quan- 
tities m and gd are mildly affected by the uncertainity of 
Z through the factor zi [see eas.([2"0" Jl -l|2"T ]l ]. 



Eas. (|34[) provides an extreme high estimate for the de- 
crease of T> and hence an extreme high estimate for Z. The 
A-body simulations show that the violent decrease of T> 



is restricted to a fa c tor of order one at each violent phase 
1 Peirani et al.1 120061; IHoffman et al.l 120071; lLapi fc Cavalierd 



120091 ; IRomano-Diaz et al.ll2006l . 120071 ; IVass et alll2009h 



In summary, the linear approximation suggests a reduc- 
tion of D at each violent phase by a factor > 1 while such 
approximation is valid. Succesive violent phases can reduce 
D by a factor u p to ~ 10 in the rang e < z < 30 as shown i n 
the simulations Peirani et al. (2006); Hoffm an et al.l (2007); 
Lapi fc Cavalierd < |2009h ; IRomano-Diaz et al. 

I (|2006l . 12007m : 

Vass et al.1 \200$ ). 



Combining the approximate decrease of T>{z) given by 
eq. (|36[) with an upper bound of a decrease by a factor ~ 100 
for the interval < z < 30 yields in the linear approximation 
the upper bound 



Z < 1.3 x 10 1 



(44) 



and we have eo. (l43[) for Z in the spherical model. The fact 
that Z in the spherical model turns to be several orders of 
magnitude below the Z value in the linear approximation 
arises from the fact that the spherical model does include 
non-linear effects and it is therefore somehow more reliable 
than the linear approximation. 

The range 1 < Z < 10000 for dSphs from A~-body simu- 
lations corresponds to realistic initial conditions in the sim- 
ulations. 

The evolutions in the two approximations considered 
(see Table I) are simple spatially isotropic expansions, fol- 
lowed by a collapse in the case of the spherical model. There 
is no possibility of non-radial motion neither of mergers in 
these approximations contrary to the case in A^-body simu- 
lations. For such reasons, the Z- values in Table I are upper 
bounds to the true values of Z in galaxies. The largest 
bound on Z yield DM particle masses below ~ 100 keV. 
Moreover, the more reliable spherical model yields 7.1 keV 
as upper bound for the DM particle mass. 

In summary, with realistic initial conditions T> will not 
decrease more than < 10000 and it is therefore fair to as- 
sume that Z < 10000 for dSphs. 



3 DM PARTICLES DECOUPLING BEING 
ULTRARELATIVISTIC 

3.1 Decoupling at Local Thermal Equilibrium 
(LTE) 

If the dark matter particles of mass m decoupled at a tem- 
perature Ti > m their freezed-out distribution function 
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only depends on 

Pc _ p P h(t) 



where 



T d T d (t) 

That is, the distribution function for dark matter particles 
that decoupled in thermal equilibrium takes the form 



f: 



equil 



P P h(t) 



T d (t) 



= f: 



equil 



where F^ qu ' 1 j s a Bose-Einstein or Fermi-Dirac distribution 
function: 



f: 



equil 



[Pc 



Notice that for eg. (|45|) in this regime: 



± 1 



(45) 



T d 



= y + O 



where y is defined by eq.([6]) and we can use as distribution 
functions 



-\equil 



(y) 



i 



e* ± 1 



(46) 



Using eqs.([SJ) and (|45p . we find then for Fermions and for 
Bosons decoupling at LTE 



9d 



3.874 eV 
2.906 eV 



Fermions 
Bosons 



(47) 



We see that for DM that decoupled at the Fermi scale: 
T d ~ 100 GeV an d g d ~ 100, m re s ults in the keV scale as al- 



ready r emarked in Bond fc Szalayl (119831 ) ; IPagels fc Primackl 
(|1982T) : iBond. Szalav fc Turner! <| 19821 ). DM particles may 
decouple earlier with T d > 100 GeV but g d is always in 
the hundreds even in grand unified theories where T d can 
reach the GUT energy scale. Therefore, eq. (|47p strongly 
suggests that the mass of the DM particles which decou- 
pled UR in LTE is in the keV scale. 

It should be not iced that the Lee- Wei nberg 
llLee fc Weinberg! Il977l: ISato fc Kobavashil 1 19771 : 
IVvsotskv. Dolgov fc Zeldovichl Il977l) lower bound as 



well as the Cowsik-McClelland l|Cowsick fc McClelland! 
19721) upper bound follow from eg . (El) as shown in 



Bovanovskv. de Vega fc Sanchez! l|2008al ) 



Computing the integrals in eq. (|13|) with the distribution 
functions eq. (|45)l yields for DM decoupling UR in LTE 



V = g 



i 

4 TV 2 

1 

8 TV 2 



C 5 (3) 



15 C J (5 ) 
C 5 (3) _ 



= 1.9625 x 10~ d 
3.6569 x 10" 3 



Fermions 
Bosons 



where C(3) = 1.2020569 ... and <(5) = 1.0369278 . 



(48) 



Inserting the distribution function eq. (|46[) into eqs. (|T9 
and (|21|) for m and g d , respectively, we obtain 



= - keV 



9d 



Z 
9 

3 I 

9 i 



0.568 
0.484 



Fermions 
Bosons 



155 
180 



Fermions 
Bosons 



(49) 



Since g = 1 — 4, for DM particle decoupling at LTE, we see 
from eq. (|49|) that g d > 100 and thus, the DM particle should 
decouple for T d > 100 GeV. Notice that 1 < zi < 10 for 
1< Z < 10000. 




Figure 1. The free-streaming wavelength today A^ s (0) in kpc 

i 

vs. the DM particle mass m in keV times gi for ultrarelativistic 
decoupling at LTE according to eq. l|51|l . 



A further estimate for the DM mass m follows by in- 
serting eq. (|48)) for V in eo. (f24)) 



t) Z 10 km/s 
9 cr s 




keV 



0.405 Fermions 
0.347 Bosons 



(50) 

Taking into account the observed values for a s and L from 
eq. (f23)l and the fact that t) < 1.6, g ~ 1 - 4, 1 < Z^ < 10, 
eg. (|50fl gives again a mass m in the keV scale as in eg. (1491) . 
Both equations (|49[) and (|50[) yield a mass larger in 17% for 
the fermion than for the boson. 

We can express the free-streaming wavelength as a func- 
tion of the DM particle mass from eqs. (|28p and (|49[) with 
the result, 



Xfs(z) = 



/keV\ 3 kpc 1_ 



V m J 



\/TT" 



7.67 Fermions 
6.19 Bosons 



We display in fig. \T\ A/ s (0) in kpc vs. m gi in keV. 



(51) 



3.2 Decoupling out of LTE 

In general, for DM decoupling out of equilibrium, the DM 
particle d istribution function takes the form 



F d (p c ) = F d 



p c m 
T~ d 'T d ' 



(52) 



Typically, thermalization is reached by the mixing of 
the particle modes and scattering between particles that 
redistributes the particles in phase space: the larger 
momentum modes are populated by a cascade whose 
front moves towards the ultraviolet akin to a direct 
cascade in turbulence, leaving in its wake a state 
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Figure 2. X+(s), X—(s), W+(s) and W-(s) as functions of s 
according to eq, l|58j l. 



of nearly LTE but with a lower temperature than 
that of equilibrium (Bovanovskv. Destri fc de Veeal 120041 : 
iDestri fc de Vegall2006r i. Hence, in the case the dark mat- 
ter particles are not yet at thermodynamical equilibrium 
at decoupling, their momentum distribution is expected to 
be peaked at smaller momenta since the ultraviolet cascade 
is not yet completed feovanovskv. Destri fc de Vegal [20041 : 
IDestri fc de Vegall2006l l . The freezed-out of equilibrium dis- 
tribution function can be then written as 



pout of LTE 



( Pc ) = F F t 



equil 



a(t) p p h{t) 



0(Pc 



Pc) 



(53) 

where £ = 1 at thermal equilibrium and £ < 1 before ther- 
modynamical equilibrium is attained. Fo ~ 1 is a normaliza- 
tion factor and p° c cuts the spectrum in the UV region not 
yet reached by the cascade. 

Inserting the out of equilibrium distribution eq. (|53p in 
the expression for the DM particle mass eq.(|5} and using 
eq. (|46[) . we obtain the generalization of eo. (|47|l for the out 
of LTE case: 

F+(°o) 



<-ld 



9 F e 



cV 



3.593 
2.695 



*•+(«) 
f-(°°) 

F_( S ) 



Fermions 
Bosons 



(54) 



where s = p®/[£ T<j] ■ Here we used eq. (|46[) and 



F±(s) 



V 2 dy 
ev ± 1 



F + (co) = - C(3) 



F_(co) =2C(3) 
(55) 



Inserting the out of equilibrium distribution eq. (|53[) into 
eqs. (|19p and (|2ip for m and gd, respectively, and using 
cq. (|46|) . we obtain the estimates 



W±(s) keV 



0.568 
0.484 



Fermions 
Bosons 



(56) 



9d - 
where 

W±(s) = 



■g* Z* e X ± (s) 



G 3 ±(s) F^(oo) 



155 
180 



Fermions 
Bosons 



(57) 



_G^(co)i^( S ) 
Here F±(s) is defined by eq. (|55p and 

G ± (s) , 



X±(s) 



G ± (s) F±(s) 



G±(oo) F±(oo) 
(58) 



, G+(co) = ^ C(5) , G_(co) = 24 C(5) 

(59) 



For small arguments s we have 
"C 5 (3) 



W + (0) = 

5* 



W-{s) = — 

25 



X+(s) a ^> 



C 3 (5) 
C 5 (3) 



27 C 3 (5) 



X-(s) 



[2025 C(3)C(5)]s 



s 1 



= 0.5732982 . 



0.4753169 . 



0.0529923 . 



0.0398856 . 



[4320 C(3) C(5)] 5 
As seen in fig. [2] 

W±(s) ^ 1 and 



X ± (s) < 1 



for 



s > 



We see from eg. (|56[) that for relics decoupling out of LTE, m 
is in the keV range. From eqs. (l56|l - (|57|l we see that both m 
and gd for relics decoupling out of LTE are smaller than if 
they would decouple at LTE. In addition, since X± (s) van- 
ishes for s — > 0, gd may be much smaller than for decoupling 
at LTE. 



We now generalize eo. (l48p for the phase-space density 
V to the out of LTE case eq.flUJ. Using eas.p^)l. and 
(|53|l we have 



V = g 



Fo 



1.9625 x 10" 



Fermions 



3.6569 x Bosons 



(60) 



Inserting now eq. (|60[l into eo. (l24p leads to the out of LTE 
generalization of the estimate for the DM particle mass m 
eq.([5T71> 



j]Z 10 



W±(s) i/^£ keV 
V Lt 



0.405 Fermions 
0.347 Bosons 

(61) 

Taking into account the observed values for a s and L from 
eq.lHSJ and the fact that 77 < 1.6, g ~ 1 - 4, 1 < < 10, 
eg. (|61 P gives again a mass m in the keV scale as in eo. (|56[) . 
Both equations (|56|l and (|6ip yield a mass larger in 17% for 



the fermion than for the boson. 



3.2.1 An instructive example: Sterile neutrinos decoupling 
out of LTE. 

We consider in this subsection a sterile neutrino v as DM 
particle decoupling out of LTE in a specific model where 
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v is a singlet Majorana fermion (g = 2) with a Ma- 
jorana mass m„ coupled with a small Yukawa-type cou- 
pling (Y ~ 10" " 



plmg 
1981; 



to a real scalar field \ (|Chikashige et al 



Gelmini fc Roncadellil 1 198 ll; ISchechter fc Valid Il982t 
Shaposhnikov fc Tkachev 2006; M cDonald fc Sahull2009i r~y- 



is more strongly coupled to the particles in the Standard 
Model plus to three right-handed neutrinos. As a result, all 
particles (except v) remain in LTE well after v decouples 
from them. 

The distribution function after decoupling of the sterile 
neutr ino v is known for small coupling Y to be (|Bovanovskvl 

mi), 



9§{v) 



where #5 (y) = ^ 



(62) 



and the coupling r is in the range 0.035 < r < 0.35 
i|Bovanovsky||2008h . 

It is interesting to compare the small (y — > 0) and large 
momenta (y — > oo) behaviour of this out of equilibrium dis- 
tribution F^(y) with the Fermi-Dirac equilibrium distribu- 
tion ea.(|4"6")l. We find 



rpequil 
r d 



(y) 



>o 2rC (|) 



CO , C 



1.341 . 



F 



equil 



(v) 



-4^0 



v^7 



Therefore, F% (y) exhibits an enhancement compared with 
the Fermi-Dirac equilibrium distribution for small (y — > 0) 
and a suppression for large momenta (y — > oo). Quali- 
tatively, the out of equilibrium distribution eq. (|53[) exhibits 
the same effect when compared to the equilibrium distribu- 
tion F^ qu ' 1 as a consequence of the incomplete UV cascade. 

We now evaluate the relevant physical quantities insert- 
ing F%(y) in the relevant equations of sec. 13.11 We find for 
m„ from eqs.flSJ and (|62p 



m v = 2.34 ^ eV , 



(63) 



which must be compared with the LTE result for fermions 
eg. ((47)) with 3 = 2. 

The phase-space density T> from eqs.{T3]) and Q62p takes 
the value, 



V = 



6 r C*(5) 
[35 7rC(7)]i 



5.627x10" r where f(7) = 1.0083493 . 

(64) 

This result is to be compared with the LTE result for 
fermions eg. (|48[) with g — 2. 

Inserting the sterile neutrino distribution function 
eq. H62p into eqs. (|19|l and (JUJ), that take into account the 
decrease of the phase-space density due to the self-gravity 
dynamics, we obtain the following mass estimates for the v 
DM particles that decoupled out of LTE, 



0.434 keV 



g d ~ ri Z* 185 . (65) 



Again, these formulas must be compared with the LTE re- 
sult for f ermions eas.(|4"9")j . g d ~ 100 corresponds to T d ~ 100 
GeV (see 
for T d in 



Kolb fc Turner) dl99Ch ) which is the expected value 
Bovanovskvl (|2008l ). 



More precisely, for the typical range 0.035 < r < 0.35, 
from eq. (|65p we find 



0.56 keV < m v Z~ * < 1.0 keV 



15<o d Z-*<. 



while for g — 2 fermions decoupling in LTE, the mass turns 
to be smaller: m = 0.48 keV and gd larger: g d Z~& — 
184 [from eq.([49]l]. 

A further estimate for m„, independent of eq. (|65[) . fol- 
lows by inserting eq. (|64[l for T> into eq. (|24p valid for a self- 
gravitating gas of DM: 



r\ Z 10 km/s 
T as 



m v ~ 0.3105 

which gives for the typical r range 
10 km/s 



keV 



(66) 



0.40 keV < m„ 



n z 



< 0.72 keV 



[Recall that 0.1 < Z~i < 1.] 



In summary the results for t he sterile neutrino decou - 
pling out of LTE in t he mod el of Chikashi ge et al] |l98lf) : 
Gelmini fc Roncadellil (Il98lf) : [Schechter fc VauT (I1982T) : 
1200^" 



Shaposhnikov fc TkachevT (2006); iBovanovskvl (2008) 
qualitatively similar to those for fermions decoupling at 
LTE. 



4 DM PARTICLES DECOUPLING BEING 
NON-RELATIVISTIC 

Particles decoupling non-relativistic at a temperature Td <C 
m are described by a freezed-out Maxwell-Boltzmann distri- 
bution function depending on 

tjt = Tf, = T d V ■ 

-id -id 



That is, 



i. 1 

^equil f 2 2 7T 2 



45 



3.1 — 1 < 

m 



25 ni (T A 
— — gaYoo — 

45 \m 

2* ivi g d Foe 



45 



3 
X 2 



(67) 



where g d is the effective number of ultrarelativistic degrees of 
freedom at decoupling, Y(t) = n(t)/s(t), n(t) is the number 
of DM particles per unit volume, s(t) their entropy per unit 
volume, x = m/Td and Yoo f ollows from the late time limit 
of th e Boltzmann equation ( Kol b fc Turner! Il990l : iBorned 
l2003h . 

For particles that decoupled NR we obtain inserting 
eo. (|67[) into the general formula for the DM particle mass 
eq.©, 



45 Q.DM pc _ 0.748 
4 7T 2 g T3 Y x ~ gY x 



eV 



(68) 



Solving the Boltzmann equa tion gives for Y x 
l|Kolb fc Turnerlll99d : lBornerll2003h . 



Y — 



45 



4^/2tt 2 9d 



(69) 
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Notice that x > 1 since the DM particles decoupled NR. Yoo 
can also be expressed in terms of cro (the thermally averaged 
total annihilation cross-section times the velocity which ap- 
pears in the Boltzmann equation) as |Kolb fc Turner]ll99d ; 
lBornedl2003h . 



Y - 1 J 45 

7T 



p d m ao Mpi 



(70) 



(We assume for simplicity S-wave annihilation). It follows 
from this relation and eq. (|68[) that 



cro = 



0. 414 10' 
~GeV r 



x 



(71) 



The trascendental equations (|68[) and (|69[) fix the values of 
m and x. They can be combined as 

x g d keV 

This equation has solutions for x > 1 provided 

m e 



(72) 



> 



i± = 0.014 



keV ' 193.5 _ 
For x = m/Td > 1 we have the analytic solution of eq. (|72p 



— = x ~ log 



193 - 5 -^)= 5 - 265 +1 °gf- vh) 

g d keV / keV/ 



We obtain the mass of the DM particle inserting the non- 
relativistic distribution function eq. (|67|) into the general for- 
mula eq. (|19[) for m 4 with the result, 

m i T i = ^^^ Z E± (73 ) 



2 7T 2 



' gd Y a 



Combining eo. (l68l) for with eq. H73[) we obtain for the 
product m Td 



VmTd = lA7 



Z\3 

— keV 

9dJ 



NR. Maxwell - Boltzmann 



(74) 

Typica l wimps are assumed to hav e m = 100 GeV and Td = 
5 GeV (Particle Data Grouu 2009). Such value for Td implies 
g d ~ 80 (|Kolb fc Turner! Il990h . Eq. d74)) thus requires for 
such heavy wimps Z ~ 10 23 well above the upper bounds 
derived in sec. [2] (see Table I). Therefore, wimps in the 100 
GeV scale are strongly disfavoured. 



We find from eos. (|13|l and (|67[l the phase-space density 
for DM decoupling NR in LTE 



V 



2n 2 
135 V3 



gd Vo, 



;.4418xl0~ 



We obtain using here the value for Yoo from eq. 

keV 



V = 0.6315 10" 



Yoo X 2 

(75) 
(76) 



Inserting this expression for T> into the general estimate for 
the DM mass eq.([24} yields, 



Vm T d ~ 0.942 



V Z\ 3 /kpc 
1d~) V~ 



10 km/s 



keV 



(77) 

As in eq. ()74p but independently from it, we reach a result 
for \JvnTd in the keV scale assuming the DM is a self- 
gravitating gas in thermal equilibrium. 



4.1 Allowed ranges for m, Td and the annihilation 
cross section oq. 

We derive here individual bounds on m, Td and <to for DM 
particles decoupling NR. 

Using that Td < m for DM particles that decoupled NR 
we obtain from eo. (l74l) a lower bound for m and an upper 
bound on Td- Furthermore, taking into account that Td > b 
eV where b > 1 or b 3> 1 for DM particles that decoupled in 
the RD era, we obtain an upper bound for m. In summary, 

v I 2 

— ' 1.47 keV < m < MeV — 
gd) b \g d 

beV <T d < (Jpj 3 1.47 keV . (78) 

Recalling that l|Kolb fc Turner! 1 19901 ) 

g d ~ 3 for 1 eV < T d < 100 keV , (79) 
and that 1 < Z < 10 4 , we see from eqs. (|78l) that 

T d < 10.2 keV . 



1.02 keV < m < -P MeV 
b 



Notice that 6 may be much larger than one with b < 
1470 (Z/g d )3 < 21960 to ensure T d < m for consistency 
in eqs. (f78")l . 



In addition, lower and upper bounds for the cross- 
section cro can be derived. From eqs. (|71[l . (|79p and x > 1 
a lower bound follows, 



ao > 0.239 10~ 9 GeV" 2 g 



(80) 



On the other hand, upper bounds for the total DM 
self-interaction cross sections ctt have been given by 
comparing X-ray, optical and lensing observations of 
th e merging of gala x y clu st ers with iV-body sim u lations 
in iMarkevitch et "all 120041'); iRandall et all d2008h: iBrada 
toPSl) (seealsolMiralda-Escudel(|2002l '); llIennawi fc Ostriker 
j2002f i; lArabadiis et al l (j2002h V: 

,2 



<7t cm 

— < 0.7 

m gr 



3200 GeV" 



Since the annihilation cross-section must be smaller than 
the total cross-section we can write the bound 



ao < 3200 m GeV" 



(81) 

Using the upper bound eq. (|78p for m yields the upper bound 
for cro 



3.32 2 
(To < ; GeV 



(82) 



This result leaves at least five orders of magnitude between 
the lower bound [eo. (|80|l ] and the upper bound for tro. The 
DM non-gravitational self-interaction is therefore negligible 
in this context. 



Exotic models where very heavy (~ 10 TeV) DM parti- 
cles are produced very late and decouple non-relativistically 
were proposed introducing two new fine tuned param- 
eters: (a) the lifetime of unstable particles (sneutrinos) 
that decay into DM (gravitinos) (b) the mass difference 
between the two particles which must be small enough 
to led to non-rel ativistic DM jCembranos et all 2005; 
IStrigari et aill2007l l It is stated in lCembranos et all (J2005) 
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and lStrigari et ah! i|2007h that such models may describe the 
observ ed phase-space density. It is stated in lBorzumati et al] 
(2008) that it is inherently difficult to fulfil all observational 
constraints in such models. 



5 CONCLUSIONS 

Our results are independent of the particle model that 
will describe the dark matter. We consider both DM parti- 
cles that decouple being NR and UR and both decoupling 
at LTE and out of LTE. Our analysis and results refer to 
the mass of the dark matter particle and the number of 
ultrarelativistic effective degrees of freedom when the DM 
particles decoupled. We do not make assumptions about the 
nature of the DM particle and we only assume that its non- 
gravitational interactions can be neglected in the present 
context (which is consistent with structure formation and 
observations) . 

In case DM particles explain the formation of galactic 
center black holes, DM particles must be fe rmions with keV- 
scale mass ijMunvaneza fc Biermann]|2006ft . 

The mass for the DM particle in the keV range is much 
larger than the temperature during the MD era, hence dark 
matter is cold (CDM). 

A possible CDM candidate in the keV scale is a 
steril e neutrino (iDodelson fc Widrowl [l99i; IShi fc Fuller! 
19991; lAbazaiian. Fuller fc Patell l200ll; lAbazaiianl l200a 



Munvaneza fc Bierm ann 2006 ; lKusenkdl2007l ) produced via 



their mixing and oscillation with an active neutrino species. 
Other pu tative CDM candidates in the keV scale are the 
gravitino dGorbunov et al . 2008; S teffenl2009l). the light neu- 
tralino l|Profumol 120081 ) and the majoron (jLattanzi fc Valid 
120071 ). 

Actually, many more extensions of the Standard Model 
of Particle Physics can be envisaged to include a DM particle 
with mass in the keV scale and weakly enough coupled to 
the Standard Model particles. 

Lyman-a forest observations provid e indirect lower 
boun ds on the masses of sterile neutrinos (|Viel et al.ll2005l . 
120071 ) while constraints from the diffuse X-ray background 
yield upper bounds on the mass of a putat i ve sterile neu- 
trino DM particle ([Dolgov fc Hansen |2002|; Watson et alj 



20061; Igovarskv et al.l 120071; iRiemer-Sorensen et all l200d 



20071 ; lLoewenstein et al.1 20091 ) . All these recent constraints 
are consistent with DM particle masses at the keV scale. 

The DAMA/LIBRA collaboration has confirmed the 
presen ce of a signal in the keV range (Bcrnabci ct al. 
l2008al ). Whether this signal is due to DM particles in 
the keV mass scale is still u nclear |Bernabei et alj l2006t 
l2008bl ; iPospelov et~al] 12008ft . On the other hand, the 
DAMA/LIBRA signals cannot be explained by a hypothet- 
ical WIMP particle with mass > 0(1) GeV since this would 
be in confl ict with previous WIMPS direct detection ex- 
Deriments llAalseth et al.1 120081; [Fairbairn & Schwctz 200' 



peri 
IHoo 



Hooper et alj|2009l ; ISavage et al.ll2008l ; I Ahmed et al.ll200 

As discussed in sec. [4] typical wimps with m = 100 GeV 
and T d = 5 GeV (|Particle Data Group|[2009l ) would require 
a huge Z ~ 10 23 , well above the upper bounds displayed in 
Table I. Hence, wimps cannot reproduce the observed galaxy 
properties. In addition, recall that Z ~ 10 23 produces from 



cq. (|28l) an extremely short A/ 3 today 

A /s (0) ~ 3.51 10" 4 pc = 72.4 AU . 

If the flyby anomaly wo uld be expla ined by DM, a keV 
scale DM mass is preferred (|Adlerll2009l ). 

Further evidence for the DM particle mass in the keV 
scale follows by contrasting the observed value of the con- 
stant surface density of galaxies to the theoretical cal- 
culation from the lineari zed Boltzmann-Vlasov equation 
jde Vega fc Sanchez! 120091 ). Independent further evidence 
for the DM particle mass in t he keV scale i s give n by 
iTikhonov etafl (|2009l ). [See also lGilmore et. all (|2007l )]. 

In summary, our analysis shows that DM particles de- 
coupling UR in LTE have a mass m in the keV scale with 
9d > 150 as shown in sec. 13.11 That is, decoupling happens 
at least at the 100 GeV scale. The values of m and gd may be 
smaller for DM decoupling UR out of LTE than for decou- 
pling UR in LTE (see sec. l3.2p . For DM particles decoupling 
NR in LTE (T d < m) we find that JmJd is in the keV 
range. This is consistent with the DM particle mass in the 
keV range. 

Notice that the present uncertainity by one order of 
magnitude of the observed values of the phase-space density 
p s /cr 3 only affects the DM particle mass through a power 
1/4 of this uncertainity according to eqs. (|19[) - (|20[ l. Namely, 
by a factor 10* ~ 1.8. 

We find that the free streaming wavelength (Jeans' 
length) is independent of the nature of the DM particle 
except for the Z factor characterizing the decrease of the 
phase-space density through self-gravity [sec. II. lj . The val- 
ues found for the Jeans' length and the Jeans' mass for m in 
the keV scale are consistent with the observed small struc- 
ture and with the masses of the galaxies, respectively. 
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